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Abstract 

Data from the PHOBOS detector have been used to study two-particle correlations over a broad 
range of pseudorapidity. A simple cluster model parameterization has been applied to inclusive 
two-particle correlations over a range of centrality for both Cu+Cu and Au+Au collisions at 
y/s NN = 200 GeV Analysis of the data for Au+Au has recently been extended to more peripheral 
collisions showing that the previously-observed rise in cluster size with decreasing system size 
eventually reaches a maximum value. Model studies have been used to quantify the significant 
effect of limited detector acceptance on the extracted cluster parameters. In the case of Au+Au, 
correlations between a trigger particle with p T > 2.5 GeV and inclusive associated particles have 
also been studied. These reveal the presence of a 'ridge' at small relative azimuthal angle which 
extends with roughly constant amplitude out to the largest relative pseudorapidity studied. The 
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large phase-space coverage of the PHOBOS detector has enabled a quantitative understanding 
of the so-called 'ZYAM' parameter used in the subtraction of the contribution of elliptic flow to 
these triggered correlations. 



Studies of single-particle distributions extracted using data from ultra-relativistic heavy ion 
collisions have provided a wealth of information concerning the global properties of these inter- 
actions. Two-particle correlations can provide more detailed information, in particular probing 
the characteristics of the particle production process itself. The PHOBOS detector at the Rela- 
tivistic Heavy Ion Collider (RHIC) has been used to study correlations in which both particles 
in a pair are chosen from the inclusive distribution 0X1 [2J as well as those in which inclusive 
particles are correlated with a trigger particle selected to have pj > 2.5 GeV Q. 

Inclusive correlations were analyzed using a simple cluster model parameterization for par- 
ticles emitted in collisions of Au+Au and Cu+Cu at -Js m = 200 GeV [2|. The average number 
of emitted particles (cluster size, K e g) and their spread in pseudorapidity (cluster width, 6) were 
extracted for interactions over a range of centrality. It was found that the cluster sizes are quite 
large and increase with decreasing system size. Comparing results for Au+Au and Cu+Cu, simi- 
lar cluster sizes are observed for collisions at the same fraction of the total inelastic cross section, 
as opposed to systems with the same number of participants. The left panel of Fig.[T]shows data 
for the cluster sizes versus centrality in Au+Au which have recently been extended to more pe- 
ripheral collisions. The rise in cluster size is seen to saturate and possibly even begin to decrease 
in the most peripheral collisions. Continuation of this trend might result in values for the cluster 
size close to that seen for p+p collisions which is shown by the gray band in the figure. 




Figure 1: (color online) Left panel: Extracted cluster sizes, corrected for acceptance effects, for Au+Au collisions at 
-^s NN = 200 GeV as a function of the fraction of the total inelastic cross section (equal to 1.0 for the most central 
collisions). Error bars include both statistical and systematic effects. The gray band shows the acceptance-corrected 
cluster size for p+p collisions. The black line is the cluster size predicted by the AMPT model. Right panel: The ratio 
of the cluster size extracted using only the particles within a limited detector acceptance over the size found using all 
particles is plotted for a variety of models as a function of the cluster width within the limited acceptance. Note that heavy 
ion collisions at y/s NN = 200 GeV typically have measured cluster widths of about 0.7-1.0. See text for discussion. 



Model calculations were analyzed to study the effect of limited detector acceptance on the 
extracted cluster parameters 0. A range of very different models was studied, including both 
simple independent cluster production and more complicated dynamical calculations. Specifi- 
cally, the analysis determined the ratio of the cluster size extracted using a limited set of particles 
over the size found using all particles. As shown in the right panel of Fig.[T] this ratio was very 
similar in all of the models when comparing systems with similar cluster width (i.e. the RMS 
of the separation in pseudorapidity between pairs of particles from a cluster). This result seems 
intuitively reasonable since for broader particle emission from a particular cluster, it is more 
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likely that one or more of the daughter particles will fall outside the detector acceptance. What 
was not expected was the relatively large magnitude of the effect even for the PHOBOS detector 
which has the largest pseudorapidity acceptance at RHIC. In order to extract the correction to 
apply to experimental data, ratios were determined as a function of the cluster width found using 
only the restricted set of particles. For the typical cluster widths seen in heavy ion collisions at 
-Js NN = 200 GeV, a detector extending over \rj\ < 3 finds only about half of K^g - 1 where K e g 
is the cluster size. Similar, although smaller, correction factors were found for the cluster width. 
Data from a detector with a pseudorapidity acceptance smaller than PHOBOS would have an 
even more limited sensitivity to the effect of this aspect of particle production. 

Correlations between a trigger particle with pj > 2.5 GeV and inclusive associated particles 
were analyzed for collisions of Au+Au at -\ls NN - 200 GeV over a range of centralities Q. A 
particularly interesting feature of preliminary triggered correlations data from other experiments 
was the so-called 'ridge', an enhancement at small relative azimuthal angle which was extended 
in relative pseudorapidity compared to the same correlations from p+p interactions Q. The 




Figure 2: (color online) Left panel: Correlated particle yield for trigger particles with pj > 2.5 GeV emitted in central 
Au+Au collisions at -^s NN = 200 GeV as a function of relative pseudorapidity. Error bars are statistical errors and 
boxes represent the systematic error due to the uncertainty in the elliptic flow subtraction. The gray band shows the 
systematic error due to the normalization between background and signal distributions. The dashed line shows the 
PYTHIA prediction for p+p at t/s = 200 GeV shifted up by a constant offset. Right panel: The normalization between 
signal and background distributions used in the elliptic flow subtraction for Au+Au collisions at ijs NN = 200 GeV in 
three centrality ranges is plotted as a function of relative pseudorapidity. See text for discussion. 



unique feature of this analysis using PHOBOS data is the very broad accessible range in relative 
pseudorapidity between the two particles. Among other results, it was discovered that this 'ridge' 
extends to the largest pseudorapidities analyzed. As shown in the left panel of Fig. [2] the data 
are consistent with a constant height of the ridge over the entire range. The dashed line in that 
panel shows the PYTHIA prediction for p+p at the same center-of-mass energy shifted up by a 
constant value of 0.25 associated particles per unit of relative pseudorapidity. 

In extracting the triggered correlation functions, it is necessary to subtract the effect of elliptic 
flow as shown in Eq. [T] where the left side is the correlated particle density. The factor B(An) 
is a normalization based on dN/dn; s(A<p, An) and b(A<p, An) are the signal (pairs from a single 
event) and background (pairs from mixed events) distributions, respectively, both normalized to 
the number of trigger events; and V{Arf) is a convolution of the elliptic flow magnitude for the 
trigger and associated particles (see Ref. [3] for details). 
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In the absence of effects which distort the relative normalization of the signal and background 
distributions, the factor a(An) would not be needed. The right panel of Fig. [2] shows the value 
of a as a function of relative pseudorapidity for Au+Au collisions in several ranges of centrality. 
The data points were found by matching the uncorrected correlation function and the predicted 
flow signal using a variety of techniques [3 1, one of which was the so-called ZYAM or zero yield 
at minimum [5 |. In previous analyses, this adjustment was treated as a somewhat arbitrary value. 
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However, the uniquely broad range of the PHOBOS detector allowed a quantitative understand- 
ing of the origin of this factor. Three features are clear in the right panel of Fig. [2] namely that 
a is typically very close to unity, that the deviation from one decreases for more central col- 
lisions, and that the dependence on relative pseudorapidity can be characterized as a Gaussian 
peak on top of a constant value. The lines in the figure are the result of a fit to the full set of data 
points (i.e. all centralities and relative pseudorapidities) using a parameterization including both 
constant and Gaussian terms. 

The component of the a factor which is independent of pseudorapidity is due to a simple 
trigger bias; events at the higher end of a centrality bin are more likely to contain the required 
trigger particle. As a result, the distribution of signal events within a given centrality bin is 
biased towards the higher multiplicity end, as compared to the background events. By sim- 
ple combinatorics, this shift results in an artificial enhancement in the number of pairs at all 
relative pseudorapidities. Both the magnitude and centrality dependence of this pseudorapidity - 
independent component of a can be reproduced quantitatively based on the widths and average 
particle multiplicities of the centrality bins used in this analysis. 

The Gaussian component is simply a reflection of the cluster-like particle-production pro- 
cesses discussed above, which impart a relative-pseudorapidity-dependent normalization be- 
tween the signal and background distributions. Both the Gaussian width and the centrality depen- 
dence of the Gaussian peak height found in fitting the a values are consistent with those found 
in the inclusive two-particle analysis. Note that the definition of the 5 parameter in Ref. |2| 
differs from the width of the Gaussian in A77 found when fitting a by a factor of \2. Com- 
paring inclusive and triggered correlations, the magnitude of the underlying cluster effect is not 
necessarily the same since in one case there is a requirement of a high transverse momentum 
trigger particle. This difference was reflected in the fit by a single overall normalization of the 
Gaussian peak height. The centrality dependence of this peak amplitude includes the trivial 
multiplicity-dependent dilution of correlated versus uncorrelated pairs. Note that the correlation 
functions used in the cluster analysis differ from those used in the triggered correlation study in 
that this multiplicity-dependent dilution is explicitly removed in the cluster analysis (see Eq. 1 
of Ref. [2 1). However, there is also the non-trivial centrality dependence of the cluster size itself 
as seen in the left panel of Fig. [T] and described in detail in Ref. Q. 

These expectations for the characteristics of the a factor were not directly included in the fit 
although they did influence the choice of functions. As one example, the Gaussian amplitude 
included an explicit (1/multiplicity) dependence. Characteristics of the data, including a varying 
centrality bin width, precluded the use of a simple analytic form, although predictions could be 
easily calculated bin by bin. The final results of the fits were compared to the predictions based 
on the effects discussed above and found to be in good quantitative agreement. As is clear in the 
right panel of Fig. [2] and as results from, among other effects, the dramatic impact of detector 
acceptance shown in the right panel of Fig. [T] this quantitative description of the various sources 
of offset between background and signal distributions in triggered correlations would not have 
been possible without the broad pseudorapidity coverage of the PHOBOS detector. 

In summary, PHOBOS data have been used to significantly enhance our knowledge of the 
nature of particle correlations in heavy ion collisions at the top RHIC energy. The data imply 
that particles are emitted from surprisingly large clusters with a non-trivial centrality dependence 
which has been recently extended towards smaller systems. The significant effect of detector 
acceptance, even in a detector as broad as PHOBOS, has been quantified. The 'ridge' at small 
relative azimuthal angle in triggered correlations has been found to extend over the full relative 
pseudorapidity range studied. Finally, the non-trivial normalization required when subtracting 
the effect of elliptic flow on two-particle correlations has been understood quantitatively. 
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